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1 Introduction 

The main subject of this paper is the role of effects due to hadronization in the 
theoretical description of inclusive r lepton decay. This paper employs the results of 
the studies of this strong interaction process elaborated over past several years [1 -6 J 
and presents some new results recently obtained in this field [7]. 

The r lepton is the only lepton which is heavy enough to decay into hadrons. This 
feature enables one to use this process in tests of Quantum Chromodynamics (QCD) 
and entire Standard Model. The theoretical description of inclusive r lepton hadronic 
decay, similarly to the case of electron-positron annihilation into hadrons, requires 
no such phenomenological models as, for example, the so-called "Parton Distribu- 
tion Functions" involved in the analysis of Deep Inelastic Scattering processes. It is 
worthwhile to note also that the experimental measurements of r lepton decay are of 
a high accuracy. But the most interesting feature of this process is that it probes the 
hadron dynamics at energies below the mass of r lepton. 

The experimentally measurable quantity here is the ratio of the total width of 
r lepton decay into hadrons to the width of its leptonic decay, which can be decom- 
posed into three parts, namely 

r(r~ —7- hadrons - u T ) 

Rt = Z s = Rt.V + Rt.A + Rt,S- (1) 

r(r _ e~ v e v T ) 

In the right-hand side of this equation, the first two terms account for the hadronic 
decay modes involving light quarks (u, d) only and associated with vector (V) and 
axial-vector (A) quark currents, respectively, whereas the last term accounts for the 
decay modes which involve strange quark. Each of the first two terms can be further 
decomposed into two parts according to the angular momentum in the hadronic rest 
frame, namely 

/V,v — K T V + K T V , K TtA — K T A + K T A . [Z) 

In what follows we shall restrict ourselves to the consideration of parts R J T ^f and R 3 T = l 
of ratio R T (□). 
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2 Theoretical description of r lepton decay 



The theoretical prediction for the quantities on hand (j2]) reads 



^T,V/A ~~ r, I Kid | ^Ewf^QCD + ^Ew) 5 (3) 



where iV c = 3 is the number of colors, |Kid| = 0.9738 ±0.0005 is Cabibbo-Kobayashi- 
Maskawa matrix element jS], S EW = 1.0194 ± 0.0050 and S' EW = 0.0010 stand for the 
electroweak corrections (see Refs. [9HTT]). and 

C = 2/! ? /f4W)4 (4) 



" '" v v VM T V ' ' Mr 

denotes the QCD contribution. Here M T = 1.777 GeV is the mass of r lepton [8], 
m v /A stands for the total mass of the lightest allowed hadronic decay mode of r lepton 
in the corresponding channel, f(x) = (1 — x) 2 (1 + 2x), and 

R v/A (s) = — lim [n v/A (s + is) - n v/A (s - is)} = - Im lim IT V/A (s + ie), (5) 

with n v/A (g 2 ) being the hadronic vacuum polarization function. In what follows the 
superscripts "V" and "A" will only be shown when relevant. 

In general, it is convenient to perform the theoretical analysis of inclusive r lepton 
decay in terms of the Adler function [12] 



"<«"> = = (•) 

In the framework of perturbation theory its ultraviolet behavior can be approximated 
by power series in the strong running coupling a s (Q 2 ) 

D{Q 2 ) * Df crt {Q 2 ) = 1 + Y! 3=1 dj [4 £) (Q 2 )] j , Q 2 -> oo, (7) 

where at the one-loop level (i.e., for I = 1) a^(Q 2 ) = 47r/(/3 In z), z = Q 2 /A 2 , 
0o = 11 — 2nf/3, A denotes the QCD scale parameter, m is the number of active 
flavors (rif = 2 will be assumed hereinafter), and d\ = 1/tt, see papers ^3HT5] and 
references therein for the details. It is worth noting also that the function R(s) (jSJ) 
and the Adler function can be expressed in terms of each other by making use of 
the following relations (see Refs. [I21CES1IIZ] for the details) 

R(s) = ± lim r E D(-C) ^ D(Q 2 ) = Q 2 Hj^^ ds . (8) 

2m e^0+ J s+i£ Q J m 2 (s + Q*y 
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In first of these equations the integration contour in the complex £-plane lies in the 
region of analyticity of the integrand. 

It is important to outline here that all the mentioned above is only valid for 
"true physical" hadronic vacuum polarization function IT phys (g 2 ) and Adler func- 
tion D phys (Q 2 ). However, as it often happens, one has to deal with their perturbative 
approximations Il pert (g 2 ) and D pert (Q 2 ), which are valid in the ultraviolet asymptotic 
only. Besides, expressions Il pert (g 2 ) and D pert (Q 2 ) are inconsistent with dispersion 
relation (jSj), which is determined by the kinematics of physical process on hand. 

Thus, one arrives at the point where the results of perturbation theory need to 
be "merged" with relevant dispersion relations. This objective can be achieved in 
the framework of "Dispersive approach" to QCD, which will be briefly overviewed 
in Sect. HJ The theoretical description of inclusive r lepton hadronic decay within 
Dispersive approach will be performed in Sect. HI whereas the analysis of this process 
within perturbative approach will be discussed in Sect. El 



3 Perturbative approach 

In this Section, we shall study the massless limit, that implies that the masses of 
all final state particles are neglected. In this case, by making use of definitions (jSJ) 
and (jHJ), integrating by parts, and employing Cauchy theorem, the quantity A QCD (jlj) 
can be represented as 

A qcd = t^ f D(M 2 e w )(l + 2e w -2e m -e m )d6, (9) 



2tt J-n 



see, e.g., papers [9| [18H20] and references therein. In general, in Eqs. floj) and (J6J) it 
is convenient to handle the leading contributions (i.e., the terms of 0-th order in the 
strong running coupling) separately from the contributions due to strong interaction, 
namely 

R( s ) = r (o)( s ) + r W( s ), D(Q 2 ) = rfW(Q 2 ) + d^(Q 2 ). (10) 

In what follows we shall restrict ourselves to the one-loop level (£ = 1). 

In fact, the only available option within perturbative approach is to directly use 
in the theoretical expression for A QCD (despite of remarks given in Sect. [2]) the per- 
turbative approximation of hadronic vacuum polarization function IT pert (g 2 ) instead 
of its unknown "true physical" expression Il phys (g 2 ). For the case of functions fl5]) 
and (j7j), this prescription eventually results in (see Ref. [7] for the details) 



r 



(0) (s) = l <— ► rf (0) (Q 2 ) = l, (11) 



Po Jo 7r(cg + 6 2 ) 
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Figure 1: Comparison of the one-loop perturbative expression A t (fl2l) (solid curves) 
with relevant experimental data (114j) (horizontal shaded bands). The leading-order 
terms of A pcrt ( fl~2l) are denoted by horizontal dashed lines. The solution for QCD 
scale parameter A (if exists) is shown by vertical dashed band. The left and right 
plots correspond to vector and axial-vector channels, respectively. 

where c = ln(M^/A 2 ) and 

A x {6) = l+2cos(#)-2cos(3#)-cos(4#), A 2 (9) = 2sin(0)-2sin(30)-sin(40). (13) 

Let us proceed now to the comparison of one-loop perturbative result ( fl2l) with 
corresponding experimental data. First of all, it is worth emphasizing here that 
perturbative approach gives identical predictions for functions A^p @ in vector 
and axial-vector channels (i.e., A pert = A pert ). However, their experimental values 
extracted from data presented in Refs. [2"UH2"2"] are different, namely 

A^ xp = 1.221 ± 0.057, A* xp = 0.748 ± 0.032. (14) 

These quantities are juxtaposed with perturbative result (1T2"]) in Fig. [TJ As one 
can infer from this figure, for vector channel the corresponding value of QCD scale 
parameter is A = (465^54) MeV (formally, there is also the second solution, A = 
(1646^29) MeV, which will not be considered hereinafter). As for the axial-vector 
channel, the perturbative approach fails to describe experimental data on r lepton 
hadronic decay, since for any value of A the function A pert f fl2|) exceeds A^ xp ffT4"|) . 

4 Dispersive approach 
4.1 General remarks 

It is crucial to emphasize that the analysis presented in Sect. |3] entirely leaves out the 
effects due to hadronization, which play an important role in the studies of strong 
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interaction processes at low energies. Specifically, the mathematical realization of 
the physical fact, that in a strong interaction process no final state hadrons can be 
produced at energies below the total mass of the lightest allowed hadronic final state, 
consists in the fact that the beginning of cut of corresponding hadronic vacuum po- 
larization function II(g 2 ) in complex g 2 -plane is located at the threshold of hadronic 
production, but not at the point q 2 = 0. Such limitations are inherently embodied 
within relevant dispersion relations, which, in turn, impose stringent physical nonper- 
turbative constraints on the quantities on hand. Obviously, these restrictions should 
certainly be accounted for when one is trying to go beyond the limits of perturbation 
theory. 

The nonperturbative constraints, which dispersion relation (jSJ) imposes on the 
Adler function ([6]), have been merged with perturbative result ([7]) in the framework 
of Dispersive approach to QCD, that has eventually led to the following integral 
representations for functions §5§ and (EJ) (see Refs. [SHUE] for the details): 

R(s) = '"M + 'ti-vlfcW*, (15) 

D(Q') = rf<°W) + ^f>rf^, (16) 



Q 2 + m 2 Jm* rK a + Q 2 a 

where 8(x) is the unit step-function (9(x) — 1 if x > and 6(x) = otherwise) 
and p(a) denotes the so-called spectral density. It is worth mentioning that in the 
massless limit (m — > 0) expressions ( Tl5|) and ( fl6l) become identical to those of the so- 
called "Analytic perturbation theory" [231426] . But it is essential to keep the hadronic 
mass m nonvanishing within the approach on hand. 

Let us proceed now to the description of inclusive r lepton hadronic decay within 
Dispersive approach. It is worthwhile to note here that there are two distinctions 
between the approach on hand and the massless perturbative approach presented 
in Sect. El Specifically, the first distinction is the incorporation of effects due to 
hadronization, and the second one is the expression for the one-loop spectral density 

A 2 

/3 hi 2 (<x/A 2 ) + vr 2 ' ~ 

which resembles previously studied nonperturbative model [T| |2"| IB"| |2"TH2"§] . In the 
right-hand side of Eq. (fTTJ) the first term is the one-loop perturbative contribution 
whereas the second term represents intrinsically nonperturbative part of the spectral 
density. 

4.2 Abrupt kinematic threshold 

There are two options in the framework of Dispersive approach. The first one is the 
so-called "abrupt kinematic threshold", which implies that the leading-order term 
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Figure 2: Comparison of expression (TT9l) (solid curves) with relevant experimental 
data ( IT4l) (horizontal shaded bands). The leading-order terms of Aq^ ( 1T91) are de- 
noted by horizontal dashed lines. The solution for QCD scale parameter A (if exists) 
is shown by vertical dashed band. The left and right plots correspond to vector and 
axial-vector channels, respectively. 



of R(s) (J5J) is approximated by the step-function 

V s J Q 2 + m^ /A 

see papers [2JIHE] and references therein for the details. Equation (ITS)) accounts only 
for basic kinematic restriction on hadronic vacuum polarization function H(q 2 ) and, 
in fact, represents a rather rough approximation. In this case the quantity Aq(,* (j4]) 
reads [6] 

where ( V/A = ml /A /M 2 and 

H{x) = g(x)0(l-x)+g(l)O(x-l)-g(( v/A ), g{x) = x(2 - 2x 2 + x 3 ). (20) 

However, similarly to perturbative case (112"]) . expression ( Tl9~l) is unable to describe the 
experimental data on r lepton hadronic decay in axial-vector channel (}Tlj) . see Fig. [2j 

4.3 Smooth kinematic threshold 

The second, more accurate, option within Dispersive approach is the so-called "smooth 
kinematic threshold". In this case the leading-order term of R(s) (JS]) takes the fol- 
lowing form [7] 

r^ A (s) = (l-^f) 3/2 d^ A (Q 2 ) = l + ^{l+ U (01n^/l + 2e[l-«(0]}, (21) 
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Figure 3: Comparison of expression (1221) (solid curves) with relevant experimental 
data (IHj) (horizontal shaded bands). The leading-order terms of Aq^ ( 1221) are de- 
noted by horizontal dashed lines. The solutions for QCD scale parameter A are 
shown by vertical dashed bands. The left and right plots correspond to vector and 
axial-vector channels, respectively. 



where = + £ 1 and £ = Q 2 /m\ /A . Eventually this leads to the following 
expression for the quantity A^, (TJJ: 



A V/A 
QCD 



1 Cv/A ( 1 + 6Cv/A n Cv/A + i R Cv/A 



1 2 1^3 



3Cv/a 1 + oCv/a ™C 



v/a 22 Sv / A 



111 



3 

16 
2 



+ 



'V/A 



p{a 



da 



a 



V/A 



(1 



1-0 



V/A 



(22) 



see paper [7] and references therein for the details. 

It is worth noting also that in the massless limit (m — > 0) both equations f fl9|) 
and (122]) acquire the same form 



A 



QCD 



a 



Ml 



da 



!+/ h[—,)p(<T)—, h(x)=g(x)9(l-x)+g(l)9(x-l). (23) 



a 



In the perturbative case the difference between expressions (I23p and (TT2"]) is due to 
the residue term 

A rcs = ^i(^), /ii(a;) = /i2(x)0(l-x) + /i 2 (l)0(x-l), h 2 (x)=x(2-2x 2 - x 3 ), 

(24) 

which appears to be additionally accounted for in Eq. (Tl2l . see discussion of this issue 
in paper [7]. 

The comparison of obtained result fl22|) with experimental data ffl4"|) yields nearly 
identical solutions for QCD scale parameter A in both channels, see Fig. |3j Namely, 
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Figure 4: The solutions for QCD scale parameter A obtained within Dispersive ap- 
proach f )22|) for vector and axial-vector channels (vertical dashed green bands) and 
perturbative solution ( [12]) corresponding to vector channel (vertical dashed light-blue 
band) . 



A = (412 ± 34) MeV for vector channel and A = (446 ± 33) MeV for axial-vector one. 
Besides, as one can infer from Fig. HJ both these solutions for QCD scale parameter 
agree very well with perturbative solution A = (465l}54) MeV obtained in Sect. [3j 

5 Conclusions 

Theoretical description of inclusive r lepton hadronic decay is performed in the frame- 
work of Dispersive approach to QCD. The significance of effects due to hadronization 
is convincingly demonstrated. The approach on hand is capable of describing exper- 
imental data on r lepton decay in vector and axial-vector channels. The vicinity 
of values of QCD scale parameter A obtained in both channels testifies to the self- 
consistency of developed approach. 
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